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Abstract 

We investigate the chiral phase transition at finite temperature and chemical 
potential within SU(2) and SU(3) Nambu-Jona-Lasinio type models. The be- 
havior of the baryon number susceptibility and the specific heat, in the vicinity 
of the critical end point, is studied. The class of the critical points is analyzed 
by calculating critical exponents. 



1 Introduction 

Strongly interacting matter at non-zero temperature and chemical potential is 
an exciting topic for physicists coming from different areas, either theoretical 
or experimental. One of the main goals in the heavy-ion physics program 
nowadays is to study the effects of several macroscopic phenomena occurring 
under extreme conditions. The discussion about the existence of a tricritical 



point (TCP) or a critical end point (CEP) is also a topic of recent interest. 
As is well known, a TCP separates the first order transition at high chemical 
potentials from the second order transition at high temperatures. If the second 
order transition is replaced by a smooth crossover, a CEP which separates 
the two lines is found. At the CEP the phase transition is of second order 
and probably falls into the same universality class of the three-dimensional 
Ising model. The existence of the CEP in QCD was suggested at the end 
of the eighties , and its properties in the context of several models have 
been studied since then HI 13 E]) ^ ^pj-^^ most recent lattice results with Nf = 
2 + 1 staggered quarks of physical masses indicate the location of the CEP at 
T^EP ^ i62±2MeV and fi^^^ = 36O±4OMeV0, however its exact location 
is not yet known. 

This point of the phase diagram is the special focus of the present contri- 
bution. Nambu-Jona-Lasinio (NJL) type models arc used and the main goal 
is to locate the critical end point and confront the results with universality 
arguments. 

We remark that most of the work done in this area has been performed 
with non strange quarks only. We will discuss the class of the critical points 
by including the analyzes in the chiral limit of both SU(2) and SU(3) versions 
of the NJL model. 

The Lagrangian of the SU(3) NJL model is given by: 

8 

C = q(^a•7-m)g + f 5][(gA''g)' + (g(z75)A>)' 



9D 



det [q{l + 75)q] + det [q{l - 75)g] . (1) 



Using a standard hadronization procedure® the baryonic thermody- 
namic potential, i^{T, V, fii), is obtained directly from the effective action. The 
baryon number susceptibility xb and the specific heat C describe, respectively, 
the response of the baryon density pB and the entropy S with respect to the 
chemical potential /Xj and the temperature T: 

X.- - E C^^C^) . (2) 



These physical quantities are relevant obscrvablcs to be studied in the 
context of possible signatures for chiral symmetry restoration. 



Our model of strong interacting matter can simulate regions of a hot and 
dense fireball created in a heavy-ion collision. Since electrons and positrons are 
not involved in the strong interaction, we impose the condition /ig = 0. So, we 
naturally get the chemical equilibrium condition /^m = yUd = /^s = /^s that will 
be used. 

After this presentation of the model, we discuss the phase diagrams in 
sec. 2. The behavior of the baryon number susceptibility and the specific 
heat in the T ~ hb plane around the CEP is studied in sec. 3, as well as the 
corresponding critical exponents. Finally, we conclude in sec. 4 with a brief 
summary of our results. 

2 Phase diagrams in SU(2) and SU(3) NJL models 

In this section we analyze the phase diagrams in different conditions in the 
T—fiB plane. Depending on the number of quark flavors Nf = 2 or Nf = 3, and 
on the masses of the quarks, different situations can occur and the transition 
from hadronic matter to QGP may be of first order, second order, or a crossover 
transition. 

We start by analyzing the differences between the three-flavor NJL model 
and its simpler version in the SU(2) sector. The phase diagrams for both models 
are presented in fig. [T]as a function of /i^ and T. 

Concerning the SU(2) model, and using physical values of the quark 
masses, rriu — nid = 6 MeV, together with A = 590 MeV and gsA^ — 2.435, 
we find that the CEP is locahzed at T<^^^ = 79.9 MeV and = 331.72 

MeV. We verified that, in the chiral limit, the transition is of second order at 
/is = and, as fiB increases, the line of second order phase transition will end 
in a first order line at the TCP. The TCP is located at = 286.1 MeV and 

rpTCP ^ 112.1 MeV. 
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Figure 1: Phase diagram in the SU(2) (left) and SU(3) (right) NJL models. 
The solid line represents the first order transition, the dashed the second order 
and the dotted the crossover transition. 

For the SU(3) NJL model, also in the chiral limit (m^ = = lUg = 0), 
we verify that the phase diagram does not exhibit a TCP: chiral symmetry 
is restored via a first order transition for all baryonic chemical potentials and 
temperatures (see right panel of fig. [T]). This pattern of chiral symmetry 
restoration remains for m„ ~ rud — and < m"'**. In our model we found 
m^"* = 18.3 MeV for to„ = to^ = 0®. When > mf^\ at fiB = 0, the 
transition is second order and, as /is increases, the second order line will end 
in a first order hue at the TCP. The TCP for = 140.7 MeV is located 
at 11%^^ = 265.9 MeV and T^^^ = 100.5 MeV. If we choose m„ = ^ 0, 
instead of second order transition we have a smooth crossover which critical line 
will end in the first order line at the CEP. Using the set of parameters ^ ^ : 
ruu = rud = 5.5 MeV, rUs = 140.7 MeV, gsA^ = 3.67, goA^ = -12.36 and 
A = 602.3 MeV, this point is localized at T'='^^ = 67.7 MeV and fi^^^ = 318.4 
MeV. 

We point out that both situations are in agreement with what is expected 
at /.is = : the phase transition in the chiral limit is of second order for 
Nf = 2 and first order for Nj >3. 

3 Susceptibilities and critical exponents in the vicinity of the CEP 



The phenomenological relevance of fluctuations around the CEP/TCP of QCD 
has been recognized by several authors. 
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Figure 2: Response functions in the SU(3) NJL model. Left panel: Baryon 
number susceptibility as a function of fj,B for different T around the CEP 
(T^EP ^ gy 7 ]^^Y ^^^y ^ rj^cEP ^ j^^Yj^ j^-gf^^ p^^^i. Specific heat 

as a function of T for different fiB around the CEP (n%^^ = 318.4 MeV and 
= ^g^^±10 MeV). 

In the left panel of fig. ^ the baryon number susceptibility is plotted for 
three different temperatures around the CEP. For temperatures below t'-^^p 
the phase transition is first order and, consequently, xb has a discontinuity. 
For T = T'-^^p the susceptibility diverges at /.jb = t^-B^^ {^^^ slope of the 
baryon number density tends to infinity). For temperatures above T'~^^^, in 
the crossover region, the discontinuity of xb disappears at the transition line. 
A similar behavior is found for the specific heat for three different chemical 
potentials around the CEP, as we can observe from the right panel of fig. O 
These calculations have been performed in the SU(3) NJL model, but the same 
qualitative behavior can be found in the SU(2) NJL version HI. 

Summarizing, the baryon number susceptibility and the specific heat di- 
verge at T = T'^^P and ^ = ^Jp^^ , respectively. 12 El El In order to make 
this statement more precise, we will focus on the values of a set of indices, the 
so-called critical exponents, which describe the behavior near the critical point 
of various quantities of interest (in our case e and a are the critical exponents 
of Xb and C, respectively). If the critical region of the CEP is small, it is 
expected that most of the fluctuations associated with the CEP will come from 
the mean field region around the CEP©. 

To a better understanding of the critical behavior of the system, we also 
analyze what happens in the SU(2) NJL model. 




Figure 3: Baryon number susceptibility as a function of — at fixed 

temperature T = T<^^^ in SU(2) (left panel) and SU(3) (right panel) NJL 
models. 

To obtain the critical exponent e(e') for the baryon number susceptibihty, 
we win consider a path parallel to the //s-axis in the T — /is plane, from lower 
(higher) /is towards the critical /i^^^, at fixed temperature T = T^^^ , To this 
purpose we consider a linear logarithmic fit of the type In^s = — In — 
Mb^^I + 5 where the term ci {c[) is independent of fiB- 

The values presented in fig. [3] for these critical exponents, calculated in 
both SU(2) and SU(3) NJL models, are consistent with the mean field theory 
prediction e = 2/3. This means that the size of the region is approximately the 
same independently of the direction of the path parallel to the ^g-axis. 

Paying now attention to the specific heat around the CEP, we have used 
a path parallel to the T-axis in the T — fiB plane from lower/higher T towards 
T'-^^p at fixed pLB = Mb^^- fig- Swe plot C as a function of T close to the 
CEP in a logarithmic scale for both SU(2) and SU(3) calculations. In this case 
we use a linear logarithmic fit, InC = —a In |T — t'-^ep\ + C2, where the term 
C2 is independent of T . 

Starting with the SU(2) case, we observe (see left panel of fig. H]), for 
T < T'^EP^ ^^^^ ^Yve slope of the fitting of data points changes for |r - rC£;P| 
around 0.3 MeV. So we have a change from the critical exponent a = 0.59±0.01 
to ai ~ 0.45 ± 0.01. As pointed out in®, this change of the exponent can be 
interpreted as a crossover of different universality classes, with the CEP being 
affected by the TCP. It seems that the effect of the hidden TCP on the CEP 
is relevant for the specific heat contrarily to what happens to xs- 
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Figure 4: Specific heat as a function of T for different values of fis around 
fj, = fJ'B^^ in SU(2) (left panel) and SU(3) (right panel) NJL models. 

We also observe that there is no clear evidence of change of the slope of 
the fitting of data points in the three- flavor NJL model (see fig. [4l right panel) . 
In fact, now we only obtain a critical exponent a = 0.61 ±0.01 when the critical 
point is approached from below. When the critical point is approached from 
above the trivial exponent a' = 0.67 ± 0.01 is obtained. 

To justify the possible effect of the hidden TCP on the CEP, as suggested 
jji El El ^ we analyze the behavior of the specific heat around the TCP. We 
find nontrivial critical exponents a = 0.40 ± 0.01 and a ~ 0.45 ± 0.01, for 
SU(2) and SU(3) cases, respectively. This result, in spite of being close, is not 
in agreement with the respective mean field value {a = 1/2). However, they 
can justify the crossing effect observed. We notice that the closest distance 
between the TCP and the CEP in the phase diagram occurs in the T-direction 
((T^TCP „ yCiJP) < (^CBP _ ^TCP))^ aj^d is ^jo^g cigar in the SU(2) case. 

4 Summary 

We verified that our model calculation reproduces qualitative features of the 
QCD phase diagram at ~ 0: for rui = the chiral transition is second-order 
for Nf = 2 and first-order for Nf > 3. Using realistic values for the current 
quark masses we find the location of the CEP in both SU(2) and SU(3) NJL 
models. 

It was shown that the baryon number susceptibility and the specific heat 
diverge at the CEP. The critical exponents for xb around the CEP, in both 



Nf = 2 and Nf = 3 NJL models, are consistent with the mean field values 
e = e' = 2/3. For the specific heat, the nontrivial values of a (1/2 < a < 2/3) 
around the CEP can be interpreted as a crossover from a mean field tricritical 
exponent (q = 1/2) to an Ising-like critical exponent (q = 2/3). 

A better insight to the difficult task of the analysis of the phase diagram of 
QCD can be provided by an extension of the NJL model where quarks interact 
with the temporal gluon field represented by the Polyakov loop dynamics. 
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